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PREFACE 


Thix  report  describes  the  software  generated  by  the  Data  Analysis  Services  group  of  RDP 
Incorporated  to  perform  the  pre*flight  and  post*flight  data  processing  and  analysis  for  the 
VIPER  mission.  Other  contributing  members  of  this  group  are  Michael  Delorey,  Peter 
Dickson,  James  Hughes,  Steven  Uicaire,  John  Palys.  Sharon  Poisson,  and  Rahul  Ruo. 

The  mission  is  l>eing  |)crformeU  by  PIAiPOB,  PI7SXA.  anil  ihcir  sup|)orting  contractors. 
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Introduction 


The  Visual  Photometric  Experiment  (VIPER)  Data  Processing  System  is 
a  series  of  programs  and  databases  that  bring  the  divergent  data 
streams  from  the  VIPER  Get-Away  Special  Canister  (GAS  Can)  to  a 
unified  form  on  a  Digital  Equipment  Corporation  (DEC)  VAX  computer 
system.  Data  streams  originating  from  a  visual  radiometer,  Xybion 
camera,  and  Pulnix  camera  in  the  GAS  Can,  NASA  shuttle  ephemeris 
tape,  and  Yale  Bright  Star  Catalogue  are  unified  by  this  software 
system  to  produce  a  spatial  database.  The  spatial  database, 
described  in  detail  later  in  this  report,  is  the  foundation  from 
which  further  zodiacal  studies  and  investigations  will  be  launched. 

The  VIPER  Data  Processing  System  can  be  divided  into  three  general 
sections.  The  first  section  consists  of  the  generation  of 
Preprocessing  databases.  These  databases  are  created  from  non-GAS 
Can  data.  The  second  section  deals  with  transporting  the  recorded 
instrument  data  from  the  GAS  Can's  two  distinct  data  recorders  to 
the  DEC/VAX  computer  system.  The  third  section  deals  with  the 
merging  of  the  separate  data  seta  from  the  first  two  sections  into 
the  spatial  database. 

•eetlen  Z  -  Preprocessing  Databases 

As  an  initial  step  before  the  Data  Processing  System  can  be  fully 
utilized,  three  databases  must  be  generated  from  data  not  gathered 
from  the  GAS  Can.  These  databases  must  be  generated  as  inputs  for 
the  final  process  described  in  Section  III  can  be  executed  but  they 
may  ba  generated  simultaneously  with  the  decommutation  processing 
in  Section  II. 

The  first  two  databases  are  derived  from  the  Yale  Bright  Star 
Catalogue  and  its  Supplement'.  The  original  Bright  Star  Catalogue 
magnetic  tape  is  run  through  a  series  of  programs  to  generate  a 
precessed  subset  of  the  catalogue.  The  subset  of  sources  is  chosen 
on  the  basis  of  the  magnitude  of  a  source,  with  the  brighter 
sources  being  selected. 

Due  to  the  spatial  resolution  of  the  VIPER  cameras,  the  instruments 
would  "see"  only  one  source  if  two  or  more  sources  were  within  20 
arcseconds  of  each  other.  Therefore,  the  subset  catalogue  was 
processed  to  combine  adjacent  stars,  and  the  resulting  "coalesced" 
stars  are  stored  in  a  FORTRAN  Direct  Access  File  (DAF) .  This 
process  has  effectively  translated  the  Yale  Bright  Star  Catalogue 
into  a  catalogue  of  stars  as  seen  by  the  VIPER  cameras. 


Hoffleit,  D.,  and  Jaschek,  C.,  (1982)  The  BRIGHT  STAR 

CATALOGUE .  4"*  Revised  Edition,  Yale  University  Observatory,  New 
Haven,  Connecticut 


I 


For  purposes  of  attitude  determination,  a  DAF  is  not  the  best  way 
to  organize  the  data.  When  attempting  to  correlate  a  catalogue 
entry  with  a  camera  image,  a  two-step  method  is  needed.  First,  a 
reference  source  in  the  camera  image  is  compared  to  a  Possible 
Match  Source  (PMS)  in  the  catalogue  by  magnitude.  Then,  the 
neighbors  of  the  reference  source  are  compared  to  the  PHS's 
neighbors.  The  corresponding  neighbors  need  to  match  in  relative- 
distance  from  the  PMS. 

It  would  be  advantageous  to  be  able  to  access  catalogue  data  by 
magnitude  and  by  nearest  neighbors.  VAX  FORTRAN  has  this  provision 
by  way  of  its  Indexed  Sequential  Access  Method  (ISAM)  files.  ISAM 
files  provide  for  the  use  of  key  fields  on  each  record  by  which  a 
particular  record  can  be  immediately  accessed  from  the  database 
regardless  of  its  position.  Furthermore,  after  a  single  keyed- 
access  acquisition  from  the  database,  the  database  is  then 
accessible  sequentially  and  the  subsequent  records  will  be  acquired 
in  ascending  order  on  that  key  field.  The  sorting  is  performed  via 
pointers  stored  in  the  file  at  the  time  of  creation. 

Two  ISAM  files  are  generated  from  the  DAF.  The  first  is  a  star- 
coalesced  ISAM  file  where  each  record  has  two  keys  by  which  it  can 
be  accessed,  ^is  database  is  called  the  coalesced  catalogue.  The 
keys  are  catalogue  source  number  and  coded  magnitude  value.  The 
catalogue  source  number  key  is  unique  to  each  source  while  the 
coded  magnit*Xxe  key  may  not  be  unique.  The  coded  magnitude  key  is 
a  value  that  includes  the  magnitude  (M)  and  uncertainty  in  the 
magnitude  (DM),  according  to  Equation  1: 

sey  m  N  X  asao  e  (  dm  x  lo  e  o.s  )  (i) 

This  key  is  then  truncated  at  the  decimal  point  to  p'-ovide  an 
integer  magnitude  value  to  be  used  as  the  coded  magnitude  key. 
Other  data  for  each  record  are  limited  to  the  celestial  coordinates 
of  right  ascension  and  declination.  These  are  epoch  1950.0  values 
precessed  to  the  date  of  the  flight  with  adjustments  for  proper 
motion.  Sources  included  in  this  database  are  selected  according 
to  a  magnitude  threshold  limit.  Only  sources  with  magnitudes  less 
than  or  equal  to  6  are  catalogued. 

The  second  ISAM  file  contains  the  neighbors  of  each  source  in  the 
coalesced  catalogue.  A  star's  neighborhood  is  defined  by  a  circle 
with  a  radius  of  13.5  degrees,  where  13.5  degrees  is  the  maximum 
field-of-view  of  the  Xybion  camera,  measured  diagonally  across  the 
video  frame.  This  neighbor  catalogue,  referred  to  as  a  proximal 
catalogue,  has  only  one  key  for  Indexed  accessed.  This  key  is 
generated  by  taking  the  reference  star's  number  (RSN)  in  the 
coalesced  catalogue,  and  the  current  neighbor's  number,  or  CNN,  and 
defining  the  key  according  to  Equation  2: 

Key  S  R8M  X  100,000  *  CNN  (2) 
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The  RSN  is  defined  so  that  the  brightest  star  has  an  RSN  of  1,  the 
next  brightest  is  2,  and  so  forth.  Multiplying  the  RSN  by  100,000 
is  done  simply  to  shift  the  digits  five  decimal  places  allowing  a 
non-overlapping  concatenation  of  the  RSN  and  the  CNN.  The  CNN, 
like  the  RSN,  is  the  ordinal  number  of  the  source  in  the  coalesced 
catalogue.  The  CNN  is  never  equal  to  the  RSN.  After  accessing  a 
record  in  the  proximal  catalogue,  the  key  can  be  decoded  to  provide 
the  CNN.  Using  the  CNN,  the  magnitude  and  celestial  coordinates  of 
the  neighbor  can  be  obtained  from  the  coalesced  catalogue.  There 
is  no  limit  set  for  the  maximum  number  of  neighbors  that  a  source 
can  have.  However,  a  neighbor  must  have  a  magnitude  less  than  or 
equal  to  the  magnitude  cutoff  used  in  the  coalesced  catalogue 
generation,  namely  6.  After  a  source  has  been  completely  processed 
for  neighbors,  a  header  record  is  written  that  uses  Equation  1  for 
its  key  with  the  same  RSN  but  the  CNN  is  set  to  zero.  Each  record 
in  the  proximal  catalogue  has  two  data  fields.  The  first  data 
field  is  the  key  mentioned  earlier.  The  value  of  the  second  field 
is  dependent  on  the  record  type.  For  header  record  types,  the 
second  field  specifies  the  number  of  neighbors  of  this  source, 
while  on  the  neighbor  record  types,  this  second  field  specifies  the 
radial  distance  from  the  reference  star  measured  in  10 '  radians. 

Figure  1  is  a  flowchart  that  indicates  graphically  the  steps 
involved  in  processing  the  Yale  Bright  star  Catalogue  to  generate 
coait>sced  and  proximal  ISAM  databases. 

The  last  preprocessing  database  that  needs  to  be  generated  is  the 
CAS  Database.  The  items  that  are  included  in  the  CAS  data  tape 
provided  by  NASA  are  selectable  by  the  agency  acquiring  the 
information.  The  actual  data  format  for  the  tape  has  not  yet  been 
obtained  from  NASA.  Therefore,  complete  descriptions  of  the 
database  and  program  that  generates  the  database  are  not  included 
in  this  report.  A  list  is  provided  in  Appendix  A  Section  III  of 
the  chosen  items  to  be  included  on  the  tape.  FORTRAN  null  modules 
for  data  acquisition  have  been  implemented  in  the  programs  where 
this  data  will  be  used  for  purposes  of  testing  other  aspects  of  the 
programs . 

Appendix  A  provides  a  more  detailed  account  of  the  data  format  of 
the  output  files  associated  with  this  processing  section.  All  the 
software  described  in  this  section  has  been  fully  tested  and  has 
been  demonstrated  to  function  properly  on  the  test  data  t?pes 
generated  by  the  VIPER  instrument. 

Seotion  II  -  night  Data  Decommutation  and  Database  Generation 

A.  From  VR8  Recorder  to  VAX  Mainframe  (Star  Cluster  File) 

Instrumentation  on  the  VIPER  GAS  Can  includes  Pulnix  and  Xybion 
video  cameras.  The  video  outputs  from  the  cameras  are  recorded  on 
a  flight-qualified,  VHS  format,  video  cassette  recorder  (VCR).  The 
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data  tape  is  a  T-160  in  length  and  recording  is  at  the  slowest 
possible  speed,  for  approximately  2  hours  and  40  minutes  of  video 
data.  At  any  time,  the  data  from  only  one  camera  are  recorded  but 
the  cameras  are  alternately  recorded  every  12  seconds.  Phillips 
Laboratory  personnel  will  remove  the  VHS  cassette  from  the  GAS  Can 
after  the  flight  for  commencement  of  the  data  playback  process. 

The  playback  hardware  system  consists  of  a  standard  VCR  whose 
output  is  fed  into  a  Data  Translation  Corporation's  0T2851  Frame 
Grabber.  The  frame  grabber  is  software-driven  from  a  Zenith  Z248 
personal  computer.  The  Z248  is  an  IBM-compatible,  80286-based 
machine  with  a  math  co-processor  and  monochrome  monitor,  utilizing 
DOS  3.1  and  Microsoft  FORTRAN  Version  5.  A  television  monitor  is 
provided  for  displaying  the  VCR  data.  Data  Translation  provided 
the  FORTRAN  library  of  routines  that  enable  user-written  FORTRAN 
programs  to  drive  the  frame  grabber  board. 

The  program  named  FRAME  has  been  written  to  acquire  a  single  frame 
of  video  data.  Keyboard  input  is  required.  The  user  is  prompted 
for  input  parameters  governing  data  display  on  the  television 
monitor,  time  track  generation  mode  (see  description  at  the  end  of 
this  paragraph),  sampling  rate  and  the  duration  of  the  program  run 
(both  are  in  real-time  seconds) .  Camera  definition  values  and 
calibration  information  are  read  from  a  file  named  SENSOR. SPC.  A 
second  input  file,  named  INDEX. HUM,  is  read  to  provide  information 
regarding  bit  placement  of  the  time  code  relative  to  the  sync 
pattern  that  precedes  it.  These  inputs  define  bit  width  (measured 
in  pixels),  confidence  limits,  number  of  bits  in  the  sync  mark, 
number  of  bits  in  the  time  code,  and  the  relative  pixel  number  of 
each  of  the  time  code  bits.  VIPER  encodes  month,  day  number,  hour, 
minute,  second,  and  hundredths  of  a  second  of  the  controller  clock, 
each  as  an  8  bit  word.  Appendix  B  provides  a  detailed  description 
of  how  these  values  are  defined.  During  VCR  data  processing,  the 
program  commands  the  frame  grabber  board  to  hold  the  current  video 
frame  in  memory.  At  this  point,  all  subsequent  video  input  to  the 
frame  grabber  board  from  the  VCR  is  ignored  until  the  acquisition 
of  another  frame  is  needed.  The  TPANE  program  then  commands  the 
frame  grabber  to  transmit  32  rows  of  pixels  to  the  PC.  Included  in 
the  first  set  of  32  rows  is  the  time  code.  The  time  code  is  a 
series  of  pixels  on  row  15  that  are  binary-coded-decimal  (BCD) 
values  of  the  time.  The  time  code  is  decoded  and  saved  for  output. 
Thresholds  for  determining  whether  the  pixel  represents  an  "on”  bit 
or  an  "off"  bit  are  camera  dependent.  Active  camera  determination 
is  performed  by  analyzing  the  peak  pixel  value  in  the  time  code. 
Because  pixel  values  in  the  tine  code  tend  to  be  higher  when  the 
Xybion  is  active,  the  threshold  for  determining  the  active  camera 
is  a  pixel  value  of  230  (on  a  scale  of  0  through  255) .  if  any  of 
the  first  6  pixels  in  the  sync  mark  exceeds  230,  the  active  camera 
is  known  to  be  the  Xybion.  The  time  code  and  active  camera 
identification  number  are  written  to  the  output  file  PIXELS. VAl.. 
If  time  track  generation  mode  has  been  selected,  the  rest  of  this 
video  frame  is  ignored  and  another  frame  is  acquired.  Otherwise, 
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processing  continues  as  described  below. 

After  time  code  determination,  the  full  set  of  pixels  that  has  been 
acquired  is  checked  for  possible  stars.  Pixel  values  are  compared 
to  a  threshold  to  determine  if  the  pixel  is  part  of  a  star  image. 
When  a  pixel  is  found  to  have  a  value  above  the  threshold,  the 
pixel's  coordinates  and  value  are  written  to  the  output  file 
PIXELS. VAL  and  the  search  continues.  After  all  32  rows  have  been 
processed,  another  set  of  32  rows  is  read  and  this  procedure 
continues  until  all  16  sets  of  32  rows  have  been  processed.  After 
completing  this  process,  an  end-of-frame  record  is  written  to  the 
output  file  indicating  that  no  more  data  exist  for  this  video 
frame . 

Program  FRAME  provides  a  significant  video  output  that  is 
selectable  at  the  start  of  the  program's  execution,  namely  the 
option  of  selecting  data  display  on  the  television  monitor.  Any 
video  frame  acquired  for  processing  is  simultaneously  displayed. 
If  the  user  selects  the  display  of  the  visual  radiometer's  (VR) 
f ield-of-view  (FOV),  A  circle  will  be  displayed  on  the  monitor 
indicating  the  extent  of  the  FOV.  Processing  of  video  data  will  be 
slightly  faster  if  data  display  is  disabled.  Figure  2  illustrates 
a  sample  data  display  on  the  television  monitor.  The  data  was 
recorded  during  a  full  VIPER  test  on  15  November  1990.  This  data 
is  of  the  Pleiades  as  seen  by  the  Xybion  camera.  The  BCD  time  code 
is  in  the  upper  left  hand  corner  in  row  14  of  the  image.  The  VR 
FOV  is  shown  as  the  oval  above  and  to  the  left  of  the  Pleiades. 
The  distortion  introduced  by  the  monitor  causes  the  circular  FOV  to 
appear  as  an  oval.  This  data  was  displayed  on  the  monitor  and 
photographed  for  this  report. 

Though  a  file  containing  pixel  values  is  now  in  existence,  the  data 
in  the  PIXELS. VAL  file  is  not  in  a  form  that  can  be  compared  to 
catalogue  sources.  Early  scenarios  for  the  VIPER  mission  proposed 
scan  rates  up  to  one  degree  per  second,  performed  as  shuttle  roll 
maneuvers  about  the  local  vertical  axis.  This  scan  rate  would  have 
produced  blurring  of  star  images  across  two  or  three  pixels,  which 
would  have  affected  the  determination  of  the  position  of  the  source 
and  also  spreads  out  the  intensity  of  the  source  over  more  than  one 
pixel.  A  follow-on  program  to  FRAME  was  written  to  read  the  pixel 
information  and  to  determine  the  star  Information  it  contains. 

Program  SOPROCSS  runs  on  the  Zenith  Z248,  described  earlier,  and 
performs  the  star  discovery  processing  (SDPROCSS) .  The  object  of 
the  star  discovery  processing  algorithm  is  to  associate  pixels 
c%«rresponding  to  the  same  object  with  each  other,  and  to  obtain 
position  and  luminosity  information  for  each  object,  for  use  in  the 
star  field  matching  algorithm  that  is  described  later.  The  pixel 
file,  PIXELS. VAL,  is  read  one  frame  at  a  time.  Clusters  of  pixels 
are  determined  that  comprise  one  source.  An  intensity-weighted 
centroid  is  determined  for  each  cluster  and  this  intensity-weighted 
centroid  is  assumed  to  be  the  source's  location  in  pixel 
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coordinates.  By  integrating  over  the  cluster  for  source 
intensities  and  using  calibration  inforroation,  the  source's 
magnitude  can  be  evaluated.  The  calibration  information  is  the 
same  as  that  found  in  the  file  SENSOR. SPC,  which  is  used  by  program 
FRAME.  After  all  the  clusters  for  a  frame  have  been  evaluated,  the 
frame's  header  record,  which  was  read  from  the  PIXELS. VAL  file,  is 
written  out  in  the  same  format,  with  the  intensity-weighted 
centroids  and  the  magnitudes  of  each  source.  These  data  are 
written  to  the  file  STARS. CLS.  This  file  is  uploaded  to  the  PL/SC 
VAX-Cluster  for  later  processing.  Figure  3  describes  this 
processing. 

SDPROCSS  can  be  transported  to  the  VAX  computer  system  with  minimal 
effort  required.  This  may  be  necessary  if  the  star  discovery 
processing  becomes  too  time  consuming  on  a  PC.  The  data  that  will 
be  transported  will  be  larger  but  the  increased  processing  speed  of 
the  VAX  may  easily  offset  this  consideration. 

The  current  mission  scenario  includes  only  relatively  slow  roll 
scans,  with  rates  of  about  8  arc-minutes  per  second.  Thus,  the 
optical  image  size  of  a  few  arc-seconds  for  point  sources  will  only 
be  marginally  affected,  and  it  may  be  possible  to  bypass  the  pixel¬ 
clustering  process.  Conversion  from  camera  intensity  to  magnitude 
will  still  be  required. 

Appendix  B  provides  file  descriptions  for  all  of  the  input  and 
output  files  used  in  the  processing  of  VCR  data  to  their  final  form 
as  a  STARS. CLS  file.  Program  FRAME  has  been  fully  tested.  Program 
SDPROCSS  has  not  been  completely  validated.  Execution  of  the 
program  provides  adequate  preliminary  results',,  but  conversion 
parameters  from  pixel  intensity  to  correspond incv  star  magnitude 
have  not  been  determined.  These  values  are  currently  being 
determined  and  will  be  incorporated  into  the  software  and  tested  by 
the  flight  of  the  VIPER  payload. 

B.  Proa  Bunstraad  Data  Recorder  to  VAX  Mainframe  (VR  Data  Pile) 

In  addition  to  the  two  video  cameras,  there  is  a  visual  radiometer, 
the  primary  VIPER  sensor.  Its  output  is  recorded  on  a  Sunstrand 
Data  Recorder  in  the  VIPER  GAS  Can.  The  high  and  low  gain  output 
voltages  of  the  visual  radiometer  are  sampled  every  one-hundredth 
of  a  second.  The  microprocessor  that  controls  the  VIPER  GAS  Can 
groups  5  sequential  12-bit  samples  from  each  gain  type  with  various 
other  GAS  Can  outputs  to  form  a  240-bit  minor  fraM  of  data. 
Twenty  minor  frames  are  grouped  together  to  form  one  major  frame  ol 
data  that  spans  one  second  of  time,  resulting  in  one  hundred  VR 
samples.  Section  Z  of  Appendix  C  illustrates  the  VIPER  Format  with 
a  description  of  each  item  in  the  major  and  minor  frames.  Major 
frames  are  written  to  the  Sunstrand  Data  Recorder's  cartridge  as 
part  of  a  larger  record.  This  cartridge  will  be  removed  by 
Phillips  Laboratory  personnel  after  the  flight. 
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At  Phillips  Laboratory,  a  DEC  PDP-11  computer  running  the  RT-11 
operating  system  is  equipped  with  a  playback  device  for  the 
Sunstrand  cartridge.  A  program  on  the  PDP-11  accesses  the  data  on 
the  cartridge,  via  direct  memory  access  techniques,  and  writes  the 
Sunstrand  data  to  an  unlabeled,  9-track,  1600  BPI  magnetic  tape. 
The  data  acquired  from  the  Sunstrand  cartridge  include  Sunstrand 
status  information.  Appendix  C,  Section  II  describes  the  format  of 
the  tape  generated  by  the  PDP-11. 

Program  REFORMAT,  on  the  VAX,  restructures  the  PDP-11  data  tape 
into  blocks  that  are  easily  transported  to  other  computer  systems. 
REFORMAT  also  verifies  the  correct  placement  of  all  sync  marks  in 
each  major  frame  (see  Appendix  C,  Section  I)  and  eliminates 
incomplete  or  erroneous  major  frames  when  producing  the  output 
file.  Appendix  C,  Section  III  details  the  revised  Sunstrand  data 
file  blocking  structure. 

From  this  revised  database  format,  the  visual  radiometer  database 
is  generated,  according  to  the  description  furnished  by  Mazzella 
and  Larson',  using  program  VRDATA,  as  shown  in  Figure  4.  Programs 
REFORMAT  and  VRDATA  have  been  fully  tested  using  sample  data. 

Section  ZZI  -  Data  Merging 

The  final  section  of  the  VIPER  Data  Processing  System  is  the 
merging  of  the  aforementioned  databases  into  a  unified  database 
with  which  analysis  of  zodiacal  data  will  be  possible. 

The  first  step  in  generating  the  unified  database  is  to  determine 
the  center  of  the  FOV  and  orientation  of  the  VIPER  cameras  and  the 
FOV  for  the  visual  radiometer.  Program  SFMATCHING  matches  the  VCR 
data  that  was  transported  to  the  VAX  in  the  file  STARS. CLS  as 
described  at  the  end  of  Section  IX,  Part  A  with  the  coalesced  stars 
and  neighbors.  The  coalesced  stars  and  neighbors  are  the  ISAM 
databases  from  Section  1  that  are  referred  to  as  the  coalesced  and 
proximal  catalogues. 

The  processing,  which  is  detailed  in  Appendix  D,  begins  with  the 
initialization  of  default  data  values  and  the  acquisition  of  user 
processing  options.  The  options  specify  the  tolerances  to  use  on 
determining  matches,  threshold  values,  and  the  maximum  scan  rate  of 
the  VIPER  GAS  Can  due  to  the  shuttle  motion.  A  frame  of  VCR  data 
is  acquired  from  the  STARS. CI3  file.  The  stars/clusters  are  sorted 
by  Increasing  magnitude.  The  brightest  star  is  chosen  to  be  the 
reference  star  and  the  relative  distances  are  calculated  for  all 
other  sources  in  the  frame.  Data  from  the  NASA  CA.S  tape,  which 
consist  of  times  and  associated  pointing  information,  are  acquired 


Mastella,  A.  and  Larson,  K.,  (30  July  1990),  Visual 
Photometric  Experiment  Data  Analysis.  Scientific  Report  Mo.  4,  OL- 
TR-90-0222,  page  9  and  Appendix  D,  AOA2319l$ 
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to  provide  an  initial  estimate  of  the  camera  pointing.  The  tape  is 
read  until  consecutive  records  are  found  with  times  that  bracket 
the  VCR  frame's  time.  Linear  interpolation  of  pointing  information 
is  performed  using  the  data  in  the  bracketing  times'  records  unless 
an  exact  match  is  found. 

A  keyed  access  is  performed  on  the  coalesced  catalogue  to  acquire 
a  star  with  a  magnitude  within  the  specified  range  of  the  reference 
star.  The  keyed  access  provides  a  possible  match  source  (W  3  -  see 
Section  I)  for  which  neighbors  are  acquired  from  the  proximal 
catalogue.  If  the  neighbors  of  the  reference  star  have  magnitudes 
and  relative  distances  that  match  the  PMS's  neighbors,  within  the 
tolerances  provided  by  the  user's  processing  options,  the  matching 
for  this  VCR  frame  is  complete.  Camera  pointing  information  and 
orientation  can  be  determined  using  the  algorithm  described  in 
Mazzella  and  Larson’.  If  the  PMS  and  its  neighbors  do  not  match, 
a  new  PMS  is  acquired  from  the  coalesced  catalogue  and  neighbors 
are  acquired  from  the  proximal  catalogue.  The  new  PMS  is  tested 
against  the  reference  star  field.  If  no  matches  can  be  found  for 
the  video  frame's  reference  star,  that  reference  is  eliminated  from 
consideration  and  the  next  brightest  source  in  the  frame  is  used  as 
the  reference  star.  Reference  elimination  continues  until  a 
termination  condition  of  either  a  match  is  found,  two  failures  to 
find  a  match,  or  no  more  references  exist  in  the  video  frame. 

When  either  of  the  termination  conditions  occur,  the  results  of 
this  frame's  processing  are  reported.  If  a  match  is  found,  the 
video  frame's  time  and  the  camera  pointing  information,  described 
in  Appendix  E,  Section  I,  are  written  out  to  the  POINTING.DAT 
output  data  file.  If  no  match  is  found,  the  video  frame's  time, 
the  brightest  star's  magnitude,  and  the  NASA  CAS  tape's  estimate  of 
the  camera  pointing  are  reported. 

Video  information  will  not  be  available  for  the  later  times  in  the 
VIPER  data  collection  interval.  When  such  is  the  case,  suitable 
dummy  values  will  be  used  as  fillers  on  the  "no  match"  output 
records.  Once  the  processing  results  for  a  video  frame  are 
reported,  subsequent  frames  are  acquired  and  similarly  processed 
until  no  more  video  frames  exist  in  the  STARS. CLS  input  file. 

The  last  step  in  the  VIPER  Data  Processing  System  involves  the 
merging  of  the  pointing  data  in  the  MATCHING.DAT  file  with  the 
visual  radiometer  database  described  in  Section  II,  Part  B.  The 
user  provides  processing  options  that  allow  for  time  sub-ranges  and 
tolerances  to  be  selected  in  addition  to  providing  visual 
radiometer  pointing  directions  in  shuttle  body  coordinates.  The 
option  of  limiting  processing  to  either  pointing  file  data  (from 


Mazzella,  A.  and  Larson  K. ,  (30  July  1990),  visual 
Photometric  Experiment  Data  Analysis.  Scientific  Report  No.  4,  GL- 
TR-90-0222,  Appendices  A  and  B,  ADA231915 


HATCHING.DAT),  or  NASA  CAS  data,  is  provided  for  evaluating  the 
usefulness  of  these  data  sets  based  on  the  reliability  and 
availability  of  the  two  pointing  sources. 

Program  SPATIAL  begins  by  acquiring  the  first  VR  data  sample  that 
falls  within  the  time  range  set  by  the  user  processing  options.  If 
pointing  data  exist  and  star  field  referencing  is  a  selected 
processing  option,  the  appropriate  data  are  acquired  from 
MATCHING.DAT  and  tested  for  validity.  Valid  pointing  data  are 
consecutive  data  records  that  bracket,  in  time,  the  VR  data 
sample's  time  such  that  the  pointing  data  times  are  within  the  time 
tolerance  of  the  VR  data  sample's  time.  Time  tolerance  is  a  user 
processing  option.  Linear  interpolation  is  performed  to  determine 
celestial  coordinates  for  the  VR  data  sample.  The  bracketing 
pointing  data's  camera  identifications  are  used  to  determine  the 
value  of  a  pointing  source  flag  included  in  the  output  described  in 
Appendix  E,  Section  II. 

If  NASA  CAS  data  exist  for  this  time  range  and  CAS  referencing  is 
a  selected  processing  option,  appropriate  data  are  acquired  from 
the  NASA  CAS  database  and  tested  for  the  same  validity  constraints 
as  the  pointing  data.  If  the  data  are  valid,  VR  pointing 
coordinates  are  cal-^ulated  from  the  CAS  data  and  the  output  flag 
identifying  the  sources  used  in  the  pointing  determination  is  set 
accordingly.  VR  pointing  coordinates  determined  from  the  star 
field  take  precedence  over  values  found  via  CAS  database  values. 
The  offset  between  CAS  and  VR  coordinates  is  expected  to  be  a 
constant  (but  initially  unknown)  transformation.  The  determination 
of  this  transformation  will  be  Isased  on  a  comparison  of  the  CAS  and 
video  camera  attitude  information  for  the  initial  period  of  the 
flight  operations,  when  both  data  sources  are  available. 

If  no  VR  pointing  coordinate  values  can  be  found,  fill-in  values 
are  placed  in  the  output  fields.  This  event  can  happen  when 
bracketing  time  values  either  do  not  exist,  fail  the  validation 
criterion,  or  are  disabled  via  processing  options.  Before  the  next 
VR  data  sample  is  processed,  results  of  the  current  sample  are 
stored  in  the  output  buffer,  to  be  written  out  as  described  in 
Appendix  E,  Section  II.  A  detailed  processing  description  for 
program  SPATIAL  is  provided  in  Appendix  F. 

Figure  5  provides  a  flowchart  description  of  the  data  merging 
process  of  the  VIPER  Data  Processing  System. 

Two  items  of  significance  are  pending  before  complete  validation  of 
the  software  can  be  achieved.  First,  VR  and  video  camera 
calibration  have  not  been  performed.  These  items  were  pending 
successful  completion  of  test  data  generation.  Second,  a  program 
to  create  a  CAS  database  of  shuttle  ephemeris  values  must  be 
created. 
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Figure  1  -  Flowchart  of  Bright  Star  Catalogue  Development 
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Figur*  2  Sample  Video  Data  Display  with  VR  Field-Of-^ lew 
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Flowchart  of  Data  Processing  Steps  on  Zenith  Z248 


Figurtt  5  -  Flowchart  of  Data  Merging  to  Form  VIPER  Spatial 
Database 
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Glossary  of  Acronyms 


BCD  - 
BPI  - 
CNN  - 

DAF  - 

DEC  - 

DOS  - 

FOV  - 
FORTRAN 

GAS  Can 

IBM  “ 
ISAM  - 

PL/GP  - 

PMS  - 
RSN  - 


VCR  - 
VHS  - 

VIPER  - 

VR  - 

VRDATA  - 


Binary  Coded  Decimal.  A  means  of  coding  a  decimal  valu^ 
using  four-bit  bytes. 

Bits  Per  Inch.  A  measurement  of  the  density  at  which 
data  are  written  on  a  tape. 

Current  Neighbor's  Number.  The  number  refers  to  the  ISAM 
catalogue  ordinal  number  associated  with  a  star.  This 
catalogue  number  references  the  current  neighbor's 
information  in  the  Bright  Star  Catalogue. 

Direct  Access  File.  A  file  organizational  scheme  used  in 
FORTRAN. 

Digital  Equipment  Corporation.  Manufacturer  of  the  VAX 
computer  system. 

Disk  Operating  System.  The  most  prevalent  operating 
system  on  IBM  personal  computers. 

Field-Of-View. 

FORmula  TRANslator.  A  programming  language  commonly  used 
for  scientific  applications. 

Get-Away  Special  Canister  used  to  carry  payloads  in  the 
shuttle  bay. 

International  Business  Machines  Corporation. 

Indexed  Sequential  Access  Method.  A  file  organizational 
scheme  used  in  VAX/VMS  FORTRAN. 

Phillips  Laboratory/Geophysics  Directorate.  The 

organizational  name  for  the  Air  Force  installation  where 
this  mission  was  initially  conceived. 

Possible  Match  Source. 

Reference  Star's  Number.  The  number  refers  to  the  ISAM 
catalogue  ordinal  number  associated  with  a  star.  This 
catalogue  number  references  this  video  frame's  reference 
star's  information  in  the  Bright  Star  Catalogue. 

Video  Cassette  Recorder. 

A  format  for  storing  analog  video  images  on  a  VCR's 
cassette  tapes. 

Visual  Photometric  ExpeRiment.  The  acronym  given  to  this 
mission. 

Visual  Radiometer.  An  photosensitive  detector  in  the 
VIPER  package.  The  data  gathered  is  in  the  visual  region 
of  the  spectrum. 

A  program  name  as  well  as  a  format  description  of  visual 
radiometer  data.  The  program  and  format  are  detailed  in 
Mazzella  and  Larson  . 
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List  of  Programs 


READ  VIPER, CAT. FOR 

A  VAX/VMS  program  to  read  the  Yale  Bright  Star  Catalogue  and 
Supplement  to  generate  a  reformatted  version  containing  only 
the  data  needed  for  the  VIPER  star  discovery  processing. 
EXTRACT-- FQl 

A  VAX/ VMS  program  to  read  the  output  from  READ_VIPER_CAT.FOR 
and  process  the  stars  according  to  their  proper  motion  to 
their  location  on  the  date  of  data  gathering.  The  stars  are 
also  tested  against  a  magnitude  threshold  so  that  only  the 
brightest  sources  are  retained. 

MKMAIWCAT.FOR 

A  VAX/VMS  program  to  read  the  output  from  EXTRACT. FOl  and 
generate  the  DAF  catalogue. 

MKCOALIDB.FOR 

A  VAX/VMS  program  that  reads  the  DAF  catalogue  from 
MKMAINCAT . FOR  and  generates  the  Coalesced  Catalogue  ISAM 
described  in  Appendix  A,  Section  I.  This  catalogue  is 
discussed  in  Section  I  of  the  report. 

MKPROXIDB.FOR 

A  VAX/VMS  program  that  generates  the  Proximal  Catalogue  ISAM 
described  in  Appendix  A,  Section  II.  This  catalogue  is 
discussed  in  Section  I  of  the  report. 

EfiAtLE..--r-Q6 

A  Zenith-PC  program  that  acquires  a  single  frame  of  video  data 
from  the  VCR  tape.  The  output  from  this  program  is  the  file 
PIXELS. VAL  described  in  Appendix  B,  Section  III. 

SDPROCSS.FOR 

A  Zenith-PC  program  that  reads  the  PIXELS. VAL  file  created  by 
FRAME. F06,  performs  the  star  discovery  processing,  and  writes 
the  file  STARS. CLS  described  in  Appendix  B,  Section  III. 
REFORMAT. FOR 

A  VAX/VMS  program  that  reads  the  VIPER  Data  Stream  from 
magnetic  tape  and  reformats  the  data  for  easier  use  on  the  VAX 
by  subsequent  programs.  See  Appendix  C  for  data  format 
specifics. 

VRDATA.F03 

A  VAX/VMS  program  that  reads  the  reformatted  VIPER  Data  Stream 
and  generates  the  VR  database. 

SFMATCHING.FOl 

A  VAX/VMS  program  that  matches  the  star  fields  recorded  on  the 
VCR  tape  and  written  (eventually)  to  file  STARS.  CLS  to 
matching  information  in  the  Yale  Bright  Star  Catalogues,  i.e. 
coalesced  and  proximal  catalogues. 

SPATIAL. FOl 

A  VAX/VMS  program  that  merges  the  pointing  information 
determined  in  the  SFMATCHING.FOl  program  with  the  VR  data  in 
the  VR  database. 
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Appendix  A 


Section  I 

Coalesced  Catalogue 
File  Structure 

The  Coalesced  Catalogue  is  an  ISAM  data  file  created  by  the 
VAX/Record  Management  System  (RMS)  utility.  It  consists  of 
fixed-length  records  of  12  bytes  (3  longwords) .  Two  fields  on 
each  record  are  defined  to  be  keys  for  Indexed  (or  keyed)  access. 
Both  keys  are  two  bytes  of  type  INTEGER.  The  first  key  is  the 
reference  star's  number  (RSN) .  The  second  key  is  the  coded 
magnitude  with  uncertainty  key  defined  according  to  Equation  1  of 
Section  I  of  the  main  report.  The  remaining  8  bytes  (2 
longwords)  are  evenly  divided  between  the  celestial  coordinate 
values  of  right  ascension  and  declination.  The  celestial 
coordinates  are  of  FORTRAN  type  REAL. 

The  following  diagram  illustrates  the  fields  on  each  data  record 
in  the  Coalesced  Catalogue. 
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Appendix  A 


Section  II 

Proximal  Catalogue 
Tile  Structure 

The  Proximal  Catalogue  is  an  ISAM  data  file  created  by  the 
VAX/RNS  Utility  with  only  one  key  for  accessing  the  data.  It  has 
fixed-length  records  of  8  bytes  (2  longwords) .  The  key  to  each 
record  is  four  bytes  of  type  INTEGER  defined  according  to 
Equation  2  of  Section  I  of  the  main  report.  The  remaining  4 
bytes  constitute  another  value  of  type  INTEGER. 

Records  are  grouped  together  to  form  a  set  of  associated 
neighbors.  For  each  entry  in  the  Coalesced  Catalogue,  there  is  a 
set  of  associated  neighbor  records  in  the  Proximal  Catalogue, 
assuming  that  the  same  magnitude  threshold  conditions  are  used 
during  both  catalogues'  generation.  A  set  contains  at  least  one 
record.  The  following  diagram  exemplifies  how  a  set  of  neighbor 
records  would  be  organized  in  the  catalogue. 
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The  number  of  neighbors  (n)  may,  in  theory,  equal  0  but  in 
practice  it  does  not  because  of  the  camera  sensitivities  and 
large  fields-of-view  and  the  cutoff  magnitude  for  the  catalogue. 
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Appendix  A 

Section  III 
CAS  Data  Items 


The  format  for  the  CAS  data  generated  by  NASA  is  selectable  by 
the  requestor  in  that  a  large  set  of  data  is  available  and  the 
requestor  can  chose  any  subset  desired.  A  minimum  set  of 
parameters  has  been  defined  for  utilization  with  the  VIPER 
pointing  determination  programs.  A  minimum  sampling  rate  of  0.5 
Hz  is  required  for  those  periods  for  which  the  visual  radiometer 
is  recording  data.  (The  standard  CAS  data  is  generated  at  1.0 
Hertz . ) 


The  following  is  the  format  and  description  of  the  values  needed 
by  the  pointing  programs. 


fifldLA? 

Fields  2-7 t 

Fields  8-13; 

Fields  14.  15: 

Fiftlda-lSr-ia: 

JLLfiAsLIl: 


Field  24; 
Fields  25-27; 
Fields  28-30: 
Fields  31-33! 


Orbit  Number 

Greenwich  Mean  Time  (Orbiter) :  year,  month, 
day,  hour,  minute,  second 

Greenwich  Mean  Time  (Mission  Control):  year, 
month,  day,  hour,  minute,  second 
Position  and  velocity  state  vectors  of 
orbiter:  Aries  true-of-date  Cartesian 
coordinates  (km;  km/ sec) 

Orbital  elements,  inclination  (degrees)  and 
right  ascension  of  ascending  node  (degrees) 
of  orbiter:  Aries  true-of-date  coordinates 
Sunrise  (l)/sunset  (0)  flag 
Attitude  flag;  0  =  telemetry;  1  = 
interpolation;  2  =  attitude  time  line 
(inertial  hold);  3  -  attitude  time  line 
(solar  inertial);  4  -  attitude  time  line 
Uocal  vertical/local  horizontal);  5  = 
attitude  time  line  (rotor) ;  6  =  default 
attitude  (Aries  1950  axes) 

Euler  angles  for  aligning  orbiter  axes  to  UVW 
axes:  yaw,  pitch,  roll  (degrees) 

Total  angular  attitude  rate;  Aries  mean  of 
1950  coordinates  (deg/sec) 

Earth  position  unit  vector;  shuttle  body  axis 
coordinates 

Position  state  vector  of  moon:  Aries  true-of- 
date  coordinates  (km) 

Position  state  vector  of  sun:  Aries  true-of- 
date  coordinates  (km) 
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Section  i 

VIPER  Camera  Specification  File 

Prograns  FRAKE  and  SDPROCSS  require  calibration  infomation  for  the 
VIPER  cameras.  This  infomation  is  provided  to  the  progress  via 
the  SENSOR. SPC  file.  Twelve  values  are  given  on  a  single  record. 
The  first  two  values  are  integer  values  for  the  intensity 
threshold.  The  intensity  threshold  is  used  to  detemine  whether  a 
pixel  is  part  of  a  star.  The  resaining  values  are  the  critical 
distance,  %fhich  defines  two  pivels  as  belonging  to  the  sase  source, 
row  and  colusn  scan  rate  cosponents,  a  saxisus  intensity  value,  the 
uncertainty  in  the  casera  intensity,  saxisus  tise  code  intensity 
used  to  identify  which  casera  is  active,  the  identification  nusber 
associated  with  the  saxisus  tise  code  intensity,  sinisus  average 
pixel  value  for  valid  tise  codes,  and  tise  code  thresholds  that 
detemine  whether  a  bit  is  "on"  for  each  casera. 

The  following  is  an  example  of  the  contents  of  a  SENSOR. SPC  file. 
The  first  line  is  the  actual  values  acquired  by  the  program.  The 
second  line  is  not  acquired  by  the  program.  It  contains  the  names 
of  the  variables  in  the  FRAME  and  SDPROCSS  programs  into  which  the 
values  are  placed.  Variables  with  more  than  one  value  are  arrays 
containing  values  for  the  Xybion  and  Pulnix  respectively. 

}  80  80j  7.  j  1.}  I.j255.}l70.  j 153  j  1  j  2000  {120  120{ 

jlnlOThrshlDCritj  A  j  B  [  Lo  !Del_L}MaxTc|MaxId|MinAve{ IBThreshj 

Pulnix  timecode  is  assumed  brighter  in  this  sample  file.  This 
configuration  was  needed  for  IB-Nov-lssO  tests. 
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Section  II 

VIPER  VCR  Time  Code 
Bit  Position  Indexing  File 

Program  FRAME  acquires  information  from  the  file  INDEX. NUM  to 
properly  decode  the  time  value  stored  in  row  15  of  the  VCR  video 
frame.  The  information  provided  in  this  file  directs  FRAME  to  the 
pixels  that  constitute  the  binary-coded-decimal  representation  of 
the  time.  The  pixel  locations  that  are  defined  are  compared  to  a 
threshold  value.  Values  above  the  threshold  are  the  binary  digit 
1.  Values  below  the  threshold  are  the  binary  digit  0. 

The  file  contains  the  following  information. 


Variable 

_ Haffig _ 

Bitwdth 

Epsilon 

Syncs kp 


NIndxNum 


IndxNum 


Description 


Bit  spacing  measured  as  pixels  per  bit 

Offset  value  used  to  fine-tune  the  time  acquisition 

routine. 

Number  of  bits  comprising  the  sync  mark  that  need 
to  be  skipped.  (SyncSkp*BitWdth) +1  is  the  first 
pixel  in  the  time  code. 

Number  of  binary  digits  in  the  time  code.  NIndxNum 
must  always  be  48  unless  modifications  are  made  to 
the  source  code  in  program  FRAME. 

Position  of  each  of  the  NIndxNum  bits.  These  are 
pixel  offset  values  from  the  end  of  the  sync  mark. 
Eight  of  these  index  values  are  on  each  of  the  six 
lines  in  the  INDEX. NUM  file. 


A  total  of  ten  records  are  in  a  standard  INDEX. NUM  file.  The 
values  that  were  found  to  produce  the  best  results  on  the  various 
laboratory  test  VCR  tapes  are  shown  in  the  following  lines.  Note 
the  structure  in  comparison  to  the  above  descriptions. 


3.287 

.0000000 

11.000 

48 


3.287 

6.574 

9.861 

13.148 

16.435 

19.722 

23.009 

26.296 

29.583 

32.870 

36.157 

39.444 

42.731 

46.018 

49.305 

52.592 

• 

55.879 

59.166 

62.453 

65.740 

69.027 

72.314 

75.601 

78.888 

82.175 

85.462 

88.749 

92.036 

95.323 

98.610 

101.897 

105.184 

108.471 

111.758 

115.045 

118.332 

121.619 

124.906 

128.193 

131.480 

134.767 

138.054 

141.341 

144 . 628 

147.915 

151.202 

154.489 

157.776 

23 


Appendix  B 


Section  111 

VIPER  Pixel  Value  and  Star  Cluster  File  Formats 

The  VIPER  PC  programs  FRAME  and  SDPROCSS  perform  data  acquisition 
from  the  VIPER  video  tape  and  determine  the  weighted  center  of  star 
clusters  on  the  video  image.  The  intermediate  file,  PIXELS. VAL, 
and  final  data  file,  STARS. CLS,  have  similar  formats  (though  the 
meaning  of  the  data  is  different).  The  file  will  contain  ''n+2” 
records  per  frame.  For  the  PIXELS. VAL  file,  the  meaning  of  "n”  is 
the  number  of  pixels  found  to  have  values  above  the  threshold.  For 
the  STARS. CLS  file,  "n”  refers  to  the  number  of  stars  seen  in  the 
video  frame. 

A  diagram  of  the  file's  contents  appears  below.  Please  note  that 
the  record  number  listed  with  the  data  below  does  not  appear  in  the 
file.  It  is  given  here  to  make  discussion  of  the  data  easier. 


Record 


Number 

j  Data  Descr lotion 

_ Example _ 

1 

j  Date,  Time,  Camera  Id 

10/21  20:30:26. 

09  1 

2 

Column  Number,  Row  Number,  Frame 

210  176 

204 

Thru 

Grabber  Intensity  (PIXELS. VAL) 

*1 

1 

Or 

2 

j  Column  Number,  Row  Number, 

210  176 

4.38 

Thru 

Magnitude  of  each  Star 

n 

,  Cluster  (for  STARS. CLS) 

n+1 

1  End-Of-Frame  Flag 

-999  -999 

-999 

APPENDIX  C 
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Section  1 
VIPER  Format 
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VIPER  Designations 


j Input 

VIPER  INPUT  FUNCTION 

!«  per 

j Input 

[Bits/  1 

j  Word  j Frame  j 

• 

|0es1g 

1  sec 

j Range 

1  Word  1 

j  #  {  N  j  Remarks 

AO 

Solar  ♦  Press  Transducer 
power  bus  1 

o-sv 

8 

16 

1 

A1 

Vis.  Rad.  (VR)  pwr  bus  Z 

1 

0-SV 

8 

17 

1 

AZ 

Hotor  pwr  bus  3 

1 

0-5V 

8 

18 

1 

A3 

Video  pwr  bus  A 

1 

0-5V 

8 

16 

2 

A4 

Shutter/Filt  Whl  CTR  pwr  bus  S 

1 

0-5V 

8 

17 

2 

AS 

Video  Tape  1ZV  bus  Aa 

1 

0-5V 

8 

18 

z 

A6 

Camera  1ZV  bus  Ab 

1 

O-SV 

8 

16 

3 

AZ 

EXP  «1SV  bus  Sa 

1 

O-SV 

8 

17 

3 

A8 

EXP  -TSV  bus  5b 

1 

O-sv 

8 

18 

3 

A9 

EXP  «5V  bus  Se 

1 

0-10V 

IZ 

A 

1 

A10 

EXP  -SV  bus  Sd 

1 

o-sv 

8 

16 

A 

All 

Shutter  Position  Status 

10 

o-sv 

8 

10 

Even 

A1Z 

Shutter  Drive  Motor 

10 

o-sv 

8 

11 

Even 

A13 

Shutter  Control 

10 

o-sv 

8 

IZ 

Even 

AH 

Shutter  Command 

10 

o-sv 

8 

10 

Odd 

ATS 

Shutter  Position  Pot 

10 

0-10V 

IZ 

3 

Even 

A16 

Filter  Position  Status 

10 

o-sv 

8 

11 

Odd 

A17 

Filter  Drive  Motor 

10 

o-sv 

8 

IZ 

Odd 

A18 

Filter  Control 

10 

o-sv 

8 

ZB 

Even 

A19 

Filter  Commend 

10 

o-sv 

8 

Z8 

Odd 

A20 

Filter  Position  Pot  Monitor 

10 

0-10V 

IZ 

A 

Even 

AZ1 

Filter  Position  Pot  Control 

10 

0-10V 

IZ 

3 

Odd 

AZZ 

Sun  Present 

1 

o-sv 

8 

17 

A 

A23 

Sun  Sensor  Analog  A 

1 

0-10V 

8 

18 

A 

AZA 

Sun  Sensor  Analog  • 

1 

0-10V 

8 

16 

S 

AZS 

Camera  i1  Iris  Position 

1 

o-sv 

8 

17 

S 

AZ8 

Caaiera  iZ  Iris  Position 

1 

o-sv 

8 

18 

5 

A27 

Caaiera  Select 

1 

o-sv 

8 

16 

6 

AZ8 

Ambient  Pressure 

1 

o-sv 

a 

17 

6 

AZ9 

Temp  Cal  Volt 

1 

o-sv 

8 

18 

6 

A30 

Temp  Sensor  T-i 

1 

o-sv 

8 

16 

7 

A31 

Temp  Sensor  T'Z 

1 

o-sv 

8 

17 

7 

A3Z 

Temp  Sensor  T-3 

1 

o-sv 

8 

18 

7 

A33 

Temp  Sensor  TH 

1 

o-sv 

8 

16 

8 

ASA 

Temp  Sensor  T'S 

1 

o-sv 

8 

17 

8 

A3S 

Temp  Sensor  T-6 

1 

o-sv 

8 

18 

8 

ASA 

Analog  Current  xIO  (CTR) 

1 

o-sv 

8 

16 

9 

A37 

Analog  Current  xlOO  (CTR) 

1 

o-sv 

8 

17 

9 

AA8 

VR  Ni  Gain 

100 

0-10V 

IZ 

1.5.9. 

13,17 

All 

SUPERCON 

AA9 

VR  lo  lUin 

100 

o-iov 

IZ 

2g6.10 
1^  111 

All 

SUPERCON 
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VIPER  Designations 


1  Input 
joesig 

VIPER  IKPUT  FUNCTION  j 

1 

#  per 

sec 

1  Input 

1 Range 

Bits/ 

Word 

word 

• 

frame  j 
#  1  Remarks 

01 

Shaft  Encoder  (S.E.)  j 

1 

1  oig 

12 

4 

5  1 

02 

Temp  1  (CTR)  lottery 

2 

Dig 

8 

16 

0,10 

03 

Teinp  2  (CTR)  VR  IPS 

2 

Dig 

8 

17 

0,10 

04 

Temp  3  (CTR)  Xyblon 

2 

Dig 

8 

18 

0,10 

05 

Te«ip  4  (CTR)  Pulnix 

2 

Dig 

8 

16 

1,11 

06 

Tcnp  5  (CTR)  Top  Plate 

2 

Dig 

8 

17 

1,11 

or 

Tettp  6  (CTR)  VMS  Recorder 

2 

Dig 

8 

18 

1,11 

o« 

Tenip  7  (CTR)  Controller  j 

2 

1  Dig 

8 

16 

2,12 

09 

Tenp  8  (CTR)  VR  Electronics 

2 

Dig 

8 

17 

2,12 

010 

Tfnip  9  (CTR)  Tape  Recorder 

2 

Dig 

8 

18 

2.12 

oil 

Tofnp  10  (CTR)  Encoder 

2 

Dig 

8 

16 

3,13 

012 

Tenp  11  (CTR)  PowerCube  Cnvtr 

2 

Dig 

8 

17 

3,13 

013 

Current  Non  (CTR) 

2 

Dig 

8 

18 

3,13 

OU 

♦5V  Hon  (CTR) 

2 

Dig 

8 

16 

4,14 

015 

♦15V  Hon  (CTR) 

2 

Dig 

8 

17 

4,14 

016 

♦28V  Hon  (CTR) 

2 

Dig 

8 

18 

4,14 

017 

■15V  Hon  (CTR) 

2 

Dig 

8 

16 

5,15 

018 

Clock  1/100  ♦  1/10  Sec  (CTR) 

2 

Dig 

8 

17 

5.15 

019 

Clock  lee  (CTR) 

2 

1  Dig 

8 

18 

5.15 

020 

Clock  Hlnutea  (CTR) 

2 

Dig 

8 

16 

6.16 

021 

Clock  Hours  (CTR)  j 

2 

j  Dig 

8 

17 

6.16 

6.16 

022 

Clock  Oay  of  Month  (CTR) 

2 

Dig 

8 

18 

023 

Clock  Month  (CTR) 

2 

Dig 

8 

16 

7,17 

024 

Tape  ttatua 

2 

Dig 

8 

17 

7,17 

025 

Cover  itatus 

2 

Dig 

8 

18 

7,17 

026 

Mode  Itatus 

2 

Dig 

16 

8,18 

029 

POi  Status  Word 

2 

Dig 

8 

17 

18 

HI 

Controller  Handshake  Monitor 

1 

1 

1 

8 

17 

16 

N2 

Tape  Roerdr  Handshake  Monitor 

5 

1 

1 

» 

17 

1.5.9, 

T ! 

NS 

PCM  Itatus  Monitor  | 

1 

1 

1 

8 

17 
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Section  11 

The  data  tape  generated  on  the  PDP-11  computer  system  contains 
alternating  status  and  data  blocks.  A  status  block  contains  3 
words  of  information.  For  purposes  of  discussion,  all  computer 
"words”  mentioned  are  assumed  to  be  16  bits  in  length.  The  first 
word  contains  the  number  of  bytes  in  the  status  block.  The  next 
word  provides  the  total  length,  in  words,  of  the  subsequent  data 
block.  The  last  word  indicates  the  number  of  bytes  of  actual  VIPER 
data  in  the  subsequent  data  block.  The  difference  between  these 
last  two  values  constitutes  the  extra  status  information  added  by 
the  tape  generating  program.  Before  the  termination  of  the  tape 
generating  program,  a  status  block  with  zero  values  for  the  last 
two  words  is  written  to  the  tape.  Because  all  data  blocks  are  the 
same  size,  all  status  blocks,  except  the  last  one,  have  the 
following  form. 


I  Word  1  }  Word  2  }  Word  3  } 

I  0006  {  4002  }  7FF0  j  (in  hexadecimal) 

The  last  status  block  appears  as  follows. 


I  Word  1  I  Word  2  }  Word  3  } 

}  0006  {  0000  I  0000  !  (in  hexadecimal) 

The  data  blocks  are  16,387  (4002  in  hexadecimal)  words  long.  The 
first  word  indicates  the  number  of  bytes  in  this  data  block.  This 
is  followed  by  32,752  (7FF0  in  hexadecimal)  bytes  of  actual  VIPER 
data  as  recorded  by  the  Sunstrand  Data  Logger.  The  data  block  ends 
with  20  bytes  of  status  information  from  the  Sunstrand  Data  Logger. 
The  data  blocks  have  the  following  form. 


1  Word  1  Word  ! 

...  j  Word  1  Word 

1 

1  •  •  • 

1  Word  j 

I  1  }  2 

1  [  7FF0  1 

7FF1  { 

}  8006  1 

.  (in 

hexadecimal) 

j  8006  !  dddd 

j  ...  1  dddd  j 

ssss  { 

• 

• 

— 

00 

o 

> 

o 

.  (in 

"dddd"  represents  a  hexadecimal  VIPER  data  item  and  "ssss" 
represents  a  hexadecimal  Sunstrand  status  item.  All  the  data 
blocks  end  with  the  hexadecimal  value  of  80A0. 


Investigation  of  VAX/FORTRAN  Input/Output  statements  determined  an 
allowable  range  of  sizes  for  records  from  magnetic  tapes  as  being 
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between  18  and  32767  bytes'.  Both  the  status  block,  whose  length 
is  6  bytes,  and  the  data  block,  whose  length  is  32774  bytes,  are 
outside  of  this  range.  Attempts  to  read  these  tape  blocks  using 
ANSI  FORTRAN  techniques  on  the  VAX  computer  system  failed. 

However,  undocumented  VAX  VMS/FORTRAN  extensions  were  discovered  by 
Real  Time  Engineering,  Incorporated.  These  extensions  were 
employed  in  a  reformatting  program  that  writes  the  data  to  a  file 
in  smaller  block  sizes  that  are  easily  transported  among  computer 
systems.  By  mounting  the  tape  as  a  "Foreign"  tape  without  header-3 
information  and  a  block  size  of  65,000,  a  FORTRAN  program  can  OPEN 
the  tape  as  an  unformatted,  variable-length-record  file.  By  adding 
a  dummy  16-bit  integer  variable  to  the  beginning  of  the  I/O  list  on 
each  READ  statement,  the  full  amount  of  data  in  each  block  can  be 
acquired  even  though  the  FORTRAN  Language  manual  states  a  limiting 
range  exists. 


VAX  FORTRAN,  Volume  2.  Language  Reference  Manual.  June  1988, 
digital  equipment  corporation,  Maynard,  Massachusetts;  for  VMS 
Version  5.0  or  higher  and  FORTRAN  Version  6.0;  page  9-9 
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Section  III 


The  VAX  data  tape  generated  by  program  REFORMAT  has  block  sizes 
that  enable  easy  data  acquisition  on  virtually  any  computer  system. 
The  tape  is  a  1600  BPI  magnetic  tape  with  fixed-length  records  of 
binary  data.  Each  record  is  18000  bytes  long  and  contains  exactly 
30  major  frames  of  VIPER  data.  A  single  record  constitutes  a 
single  tape  block.  If  program  REFORMAT  finds  anything  amiss  in  the 
data  (i.e.  missing  sync  marks  or  incomplete  major  frames),  600 
bytes  of  zeroes  (one  major  frame)  are  written  to  the  output  stream 
and  the  corrupted  input  data  frame  is  skipped. 
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Appendix  D  -  VIPER  star  Field  Matching  Processing  Description 

The  following  is  a  detailed  description  of  the  processing  steps 
of  the  Star  Field  Matching  program.  The  pseudo-code  given  here 
is  the  basis  of  the  SFMATCHING  program.  Most  of  the  items  in 
parenthesis  are  either  the  actual  variable  names  used  in  the 
program,  names  of  files  that  are  accessed,  or  are  FORTRAN 
variable  types  associated  with  items,  items  that  are  enclosed  in 
curly-brackets  (i.e.  {))  are  module  names.  Module  names  that  are 
underlined  denote  the  beginning  of  the  pseudo-code  associated 
with  the  name.  A  module  name  imbedded  in  a  pseudo-code  statement 
is  treated  as  a  FORTRAN  call  to  that  module. 


Xiqpttt: 

1)  Star  field  list  from  PC  star  discovery  processing,  with  stars 
sorted  by  increasing  magnitude:  (Filename:  STARS. CLS) 

Frame  time.  Camera  ID,  Number  of  stars  (NumFldStar) ; 

Pixel  row  (lYPixel),  Pixel  column  (JXPixel),  and  Magnitude 
(FieldMag)  for  each  star. 

2)  Star  catalogue 

Binary  search,  by  magnitude,  is  performed  on  Coalesced  Star 
Catalogue  (COAL.IDB);  neighboring  stars  are  acquired  from 
Proximal  Catalogue  (PROX.ZDB). 

Note:  Coalesced  Star  Catalogue  is  assumed  to  be  acquired 
completely  into  memory  in  magnitude-order. 

3)  CAS  data  estimate  (see  Appendix  A  section  II  for  details) . 

frooess: 

< ACQUIRE!  Acquire  comparative  pointing  positions: 

a)  CAS  tape  result,  for  corresponding  time: 

Is  there  a  CAS  tape  result  for  this  frame  time  (or  an  adjacent 
bracketing  pair  for  interpolation)? 

Yes:  Set  condition  to  check  index  star  versus  CAS; 

No:  Clear  condition  to  check  index  star  versus  CAS; 

b)  Previous  star  field  result,  if  time  difference  is  within 
limits  (OeltaTFrame  <  DeltaTime„J  ; 

ISfiTFRAHE.)  Set  up  star  field  frame: 

Reference  "plate  scale"  conversions,  using  coordinates 

referenced  to  center  of  field,  based  on  camera  identifier: 
For  i  »  1  to  NumFldStar 

DistX(i)  =  FactX  *  (XMidPixel  -  JXPixel (i)) 

DistY(i)  =  FactY  *  (YMidPixel  -  lYPixel(i)) 

Next  i 

Define  brightest  star  magnitude,  FieldMagO: 

Initialize  error  flag  for  catalogue  match  failures  for  index 
star: 
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IndexError  =  0 
Initialization  test: 

Is  time  difference  between  current  frame  and  previous  frame 
(DeltaTFrame)  less  than  (frame  length) / (estimated  scan 
rate) ? 

["Estimated  scan  rate"  is  nominally  1  degree  per  second, 
and  will  be  defined  as  a  program  variable.] 

Yes:  Is  current  FieldMagO  approximately  equal  to  previous 
FieldMagO? 

Yes:  Use  proximal  catalogue  field  for  previous  frame  as 
initial  reference  for  current  frame;  proceed  to 
calculation  of  relative  field  star  distances 
(RELDIST); 

No:  Prepare  to  search  index  for  new  match  to  FieldMagO 
(FIRSTSTAR); 

No:  Prepare  to  search  index  for  new  match  to  FieldMagO 
(FIRSTSTAR); 

( FIRSTSTAR )  Find  index  star  with  magnitude  (CatNag)  within 
tolerances  of  FieldMagO; 

Binary  search  for  star  in  index  catalogue: 

IndexMin  0 

IndexMax  -  Number  of  index  catalogue  entries 
Set  initial  distance  comparison  counter  limit: 

iniDist  «  2 

iSSMSai 

Index  »  (IndexMin  IndexMax) /2 
If  CatMag( Index)  «  FieldMagO  Then 

Proceed  to  find  first  star  in  index  catalogue  with 
this  magnitude  (FINDFIRST); 

Else  If  CatMag( Index)  >  FieldMagO  Then 
IndexMax  »  Index 
Else 

IndexMin  «  Index 
Endif 

If  IndexMax  >  IndexMin  Then 

Continue  search  process  (SEARCH) 

Else 

[The  field  star  magnitude  occurs  in  a  gap  in  the 
catalogue  magnitudes]; 

Checic  the  catalogue  entry  found  against  the  field 
star,  with  error  tolerances; 

If  I CatMag( IndexMin)  -  FieldMagO |  <  TolMag  Then 
Set  Index  =  IndexMin 

Proceed  to  search  for  first  index  catalogue  entry 
with  magnitude  equal  to  CatMag( Index) 
(FINDFIRST) ; 

Else  If  |CatMag( IndexMin+1)  -  FieldMagO]  <  TolMag 
Then 
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Set  Index  »  IndexMin  +  1 

Proceed  to  search  for  first  index  catalogue  entry 
with  magnitude  equal  to  CatMag( Index) 
(FINDFIRST);  (it  should  be  this  entry,  the 
first  catalogue  magnitude  above  the  gap) ; 

Else 

Warning:  Error  finding  index  star! 

Set  FleldMagO  »  FieldMag(2)  (skip  to  next 
brightest  star  in  field) ; 

Make  the  original  index  star  "disappear": 
FieldMag(l)  «  FieldMag(2) 

DistX(l)  »  DistX(2) 

Disty(l)  DistY(2) 

Set  initial  distance  comparison  counter  limit: 
IniDlst  «  3 

Increment  error  counter: 

IndexError  «•  IndexError  -4-  1 
If  IndexError  <  2  Then 

Proceed  to  search  for  new  reference  star  in 
index  catalogue  (FIRSTSTAR); 

Else 

Report  failure  of  frame  match  (FAILURE); 
Proceed  to  next  star  field  (ACQUIRE); 

Endif 

Endif 

Endif 

t FINDFIRST ;  Find  first  star  in  catalogue  with  this  magnitude: 

If  Index  >  1  Then 

If  CatMag(Index-l)  »  CatNag( Index)  Then 
Index  »  Index  ->  1 

Loop  to  check  previous  entry  in  catalogue  (FINDFIRST) 
Endif 
Endif 
Set: 

IndexLow  »  index 
IndexHigh  »  Index 
Advance  »  TRUE 
RefMag  -  CatMag( Index) 

CheckLimit  =  0 

(for  use  in  finding  'next*  index  candidate,  if  necessary); 

!RELDIST>  Define  distances  from  brightest  star; 

XO  -  DistX(l) 

YO  =  DistY(l) 

For  i  »  IniDist  to  NumFldStar 

FieldDist(i)  =  sqrt  ((DistX(i)  -  XO)'  +  (DistY(i)  -  YO)’) 
Next  i 
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(PATHCHECKt  Preli ainary  check  of  brightest  field  star  against  CAS 
position; 

If  CAS  data  is  available  Then 

Compare  star  position  to  CAS  data; 

If  discrepancy  exceeds  tr'erances  Then 
Search  for  new  index  ..car  {NEXTINDEX) 

Endif 

Endif 

(CHECKFRAME)  Conpare  proximal  catalogue  entries  to  field  stars: 
Select  subset  from  proximal  catalogue  for  match  to  field 
stars: 

Clear  ProxUsed  array  (TRUE/FALSE  flags,  indicating  that  a 
proximal  catalogue  star  has  already  been  matched  to  a 
field  star):  ProxUsed(j)  =  FALSE,  for  all  j; 

BadMatch  =  0  (counter  for  number  of  unmatched  field  stars) 
For  i  =  IniOist  to  NumFldStar 
For  j  =  1  to  NumProx 

If  ProxUsed(j)  =  TRUE  Then  Next  j 
If  {FieldDist(i)-ProxDist( j) j  <  TolDist  Then 
I f  ! FieldMag ( i ) -ProxMag ( j ) j  <  TolMag 
IndexProx(i)  =  j 
ProxUsed (■})  -  TRUE 
Next  i 
Endif 
Endif 
Next  j 

Warning:  Matching  star  not  found!  [Have  completed  'j' 
loop  with  no  match.] 

If  i  <  i„.  Then 

[Star  field  shows  early  discrepancy  from  proximal 
catalogue] ; 

Search  for  new  index  star  (NEXTINDEX) 

Else 

BadMatch  =  BadMatch  +  1 

IndexProx(i)  =  0  (flag  field  star  as  having  no 
catalogue  match) ; 

If  BadMatch  >  BadMatch,„  Then  search  for  new  index 
star  (NEXTINDEX) 

Endif 
Next  i 

( POINTING )  Determine  field  center,  orientation,  and  visual 
radiometer  pointing; 

Determine  coordinate  transformation  from  camera  to  celestial 
coordinates  using  least-square  minimization; 

Can  determine  orientation  about  line-of-sight  from  one  of  the 
Euler  angles  for  the  transformation; 

Use  transformation  matrix  to  compute; 
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Celestial  coordinates  for  center  of  camera  field; 

Celestial  coordinates  for  center  of  visual  radiometer, 
based  on  its  camera  coordinate  offset; 

(See  Mazzella  and  Larson*  for  separate  description  for 
transformation  determination.) 

If  OeltaTFrame  <  DeltaTime^  Then 

Compute  angular  distance  from  previous  position  (DeltaAng) ; 
If  DeltaAng/DeltaTFrame  >  Critical  scan  rate  Then 
Warning:  Excessive  scan  rate! 

Report  results  (REPORT); 

Endif 

Endif 

(REPORT!  Report  results: 

Report  (to  pointing  file) : 

Frame  time; 

Camera  pointing  (center)  in  celestial  coordinates; 

Camera  orientation,  as  azimuthal  degrees  from  North  about 
pointing  direction; 

Root  mean  sguared  distance  error  for  star  field  versus 
catalogue; 

Root  mean  sguared  intensity  error  for  star  field  versus 
catalogue; 

Visual  radiometer  pointing  in  celestial  coordinates; 

Number  of  stars  in  visual  radiometer  field; 

Estimated  total  intensity  of  stars  in  visual  radiometer 
field; 

( SAVEFRAME )  Save  parameters  of  current  frame  for  reference  by 
subsequent  frame  (only  for  successful  field  matches) ; 

Frame  time; 

FieldMagO; 

Camera  pointing  position  in  celestial  coordinates; 

Proceed  to  next  star  field  (ACQUIRE); 

(NEXTINDEX)  Procedure  to  select  the  next  candidate  from  the  index 
catalogue  for  the  field  reference  star; 

If  IndexHigh  <  number  of  stars  in  index  catalogue  Then 
If  CatMag(Index+l)  =  RefMag  Then 
IndexHigh  =  IndexHigh  +  1 
Index  =  Index  +  1 

Proceed  to  preliminary  check  for  proximal  catalogue 
match  { PATHCHECK ) ; 

Endif 

Else  If  Advance  =  TRUE  Then 
Advance  =  FALSE 

If  IndexHigh  =  Number  of  index  catalogue  entries  Then 
CheckLimit  =  CheckLimit  +  1 


38 


Appendix  D  -  VIPER  Star  Field  Matching  Processing  Description 

Loop  on  search  for  next  index  star  { NEXTINDEX  > ; 

Endif 

If  |CatMag(IndexHigh+l)  -  FieldMagO|  <  TolMag  Then 
IndexHigh  ==  IndexHigh  +  1 
Index  =  IndexHigh 
CheckLimit  =  0 

Proceed  to  preliminary  check  for  proximal  catalogue 
match  (PATHCHECK); 

Else 

CheckLimit  =  CheckLimit  +  1 
Endif 
Else 

Advance  =  TRUE 
If  IndexLow  =  1  Then 

CheckLimit  =  CheckLimit  +  1 

Loop  on  search  for  next  index  star  (NEXTINDEX); 

Endif 

If  jCatMag(IndexLow-l)  -  FieldMagOj  <  TolMag  Then 
IndexLow  =  IndexLow  -  1 
Index  =  IndexLow 
CheckLimit  =  0 

Proceed  to  preliminary  check  for  proximal  catalogue 
match  (PATHCHECK); 

Else 

CheckLimit  =  CheckLimit  +  x 
Endif 
Endif 

If  CheckLimit  <  2  Then 

(at  least  one  end  of  search  is  still  within  magnitude 
tolerances) 

Loop  on  search  for  next  index  star  (NEXTINDEX); 

Else 

Warning;  No  further  index  star  candidates! 

Set  FieldMagO  =  FieldMag(2)  (skip  to  next  brightest  star  in 
field) ; 

Make  the  original  index  star  "disappear”; 

FieldMag(l)  =  FieldMag(2) 

DistX(l)  =  DistX(2) 

DistY(l)  =  Disty(2) 

Set  initial  distance  comparison  counter  limit; 

IniDist  =  3 

Increment  error  counter: 

IndexError  =  IndexError  +  1 
If  IndexError  <  2  Then 

Proceed  to  search  for  new  reference  star  in  index 
catalogue  ( FIRSTSTAR) ; 

Else 

Report  failure  of  frame  match  (FAILURE); 

Proceed  to  next  star  field  (ACQUIRE); 

Endif 
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End  if 

(FAILURE )  Report  failure  of  star  field  match; 

Frame  time; 

CAS  tape  estimate  for  pointing,  if  available; 

Magnitude  of  brightest  field  star; 
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Appendix  E  -  File  Format  for  Data  Merging  Programs 


Section  I  **  Output  from  the  SFNATCHING  program 

The  star  field  matching  program,  SFMATCHING,  determines  the 
attitude  of  the  VIPER  video  cameras,  and  thereby  the  visual 
radiometer,  by  comparing  a  vcr  video  frame  with  known  sources  from 
the  Yale  Bright  Star  Catalogue.  When  a  successful  match  occurs 
between  the  video  frame  and  the  catalogue,  the  following 
information  is  written  to  the  ASCII  file  P0INTING.DAT. 

A.  Time  stamp  of  frame  -  Julian  date,  hours,  minutes,  and 
seconds 

B.  Camera  Identification  number 

C.  Estimated  total  intensity  of  stars  in  VR  field 

D.  Camera  pointing  -  Celestial  coordinates  of  the  center  ol 
the  FOV  of  the  video  camera. 

E.  VR  pointing  >  Celestial  coordinates  of  the  center  of  the 
FOV  of  the  visual  radiometer. 

F.  Camera  Orientation  -  Euler  rotation  angles  from  an 
initially  aligned  camera/celestial  coordinate  system. 
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Section  IX  -  Output  from  the  SPATIAL  program 


The  spatial  data  base  generating  program,  SPATIAL,  produces  a  data 
file  that  merges  shuttle  ephemeris  data  with  the  VR  data  values. 
The  data  are  written  out  in  variable  length  blocks  with  a  maximum 
size  of  3242  bytes.  The  first  word  (i.e.  2  bytes)  contains  the 
number  of  records  in  the  block.  Each  record  is  54  bytes  in  length 
and  a  maximum  of  60  records  is  allowed  for  each  block. 


A. 

B. 

C. 

D. 

E. 

F. 

G. 

H. 

I. 

J. 

K. 

L. 

M. 

N. 

O. 
P. 


VIPER  day  count  since  1  January  1985. 

Time  of  day  in  seconds. 

Shutter  position. 

Filter  position. 

VR  temperature,  in  counts. 

VR  radiance,  in  counts. 

VR  right  ascension,  in  degrees. 

VR  declination,  in  degrees. 

Shuttle  right  ascension,  in  degrees. 
Shuttle  declination,  in  degrees. 

Shuttle  altitude,  in  km. 

Solar  right  ascension,  in  degrees. 

Solar  declination,  in  degrees. 

Lunar  right  ascension,  in  degrees. 

Lunar  declination,  in  degrees. 

Source  flag. 

0  -  no  position 

1  -  Xybion  only 

2  -  Pulnix  only 

3  >  Xybion  and  Pulnix  without  CAS 

4  -  CAS  only 

5  -  Xybion  and  CAS 

6  Pulnix  and  CAS 

7  >  Xybion,  Pulnix,  and  CAS 

Note:  Absence  of  CAS  data  implies  missing 
and  moon,  and  missing  orbital  parameters. 


4  bytes  of  INTEGER 
4  bytes  of  REAL 
2  bytes  of  INTEGER 
2  bytes  of  INTEGER 
2  bytes  of  INTEGER 
2  bytes  of  INTEGER 
4  bytes  of  REAL 
4  bytes  of  REAL 
4  bytes  of  REAL 
4  bytes  of  REAL 
4  bytes  of  REAL 
4  bytes  of  REAL 
4  bytes  of  REAL 
4  bytes  of  REAL 
4  bytes  of  REAL 
2  bytes  of  INTEGER 


coordinates  of  sun 
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The  following  is  a  detailed  description  of  the  processing  steps 
of  the  Spatial  Data  Base  Generating  program.  The  pseudo-code 
given  here  is  the  basis  of  the  SPATIAL  program.  Most  of  the 
items  in  parenthesis  are  either  the  actual  variable  names  used  in 
the  program,  names  of  files  that  are  accessed,  or  are  FORTRAN 
variable  types  associated  with  items.  Items  that  are  enclosed  in 
curly-brackets  (i.e.  ())  are  module  names.  Module  names  that  are 
underlined  denote  the  beginning  of  the  pseudo-code  associated 
with  the  name.  A  module  name  imbedded  in  a  pseudo-code  statement 
is  treated  as  a  FORTRAN  call  to  that  module. 


Inputs: 

1)  VRDATA  ■  Visual  Radiometer  data  base,  in  VRDATA  format; 

2)  NATCHING.DAT  "  pointing  data  from  star  field  matching; 

3)  CAS  data  base  ■  pointing  transformation  angles,  derived  from 
CAS  data,  with  additional  celestial  parameters. 

User  Psooessisg  Options: 

1)  Starting  month/day/year/ time  for  processing; 

2)  Ending  month/day/year/ time  for  processing; 

3)  Pointing  directions  (zenith:  VRZEN,  azimuth:  VRAZ  angles) 
for  Visual  Radiometer  in  shuttle  body  coordinates;  (VRZEN  ■> 
200  if  no  estimate  for  usage  with  CAS  is  available;  default 
is  VRZEN  «  180,  VRAZ  ■  0]  (For  shuttle  body,  X  is  forward,  y 
is  over  right  wing,  and  Z  is  toward  bottom. ) 

4)  Time  tolerance  (TIMEINT) ,  in  seconds,  for  bracketing 
pointing  information  relative  to  visual  radiometer  data  time 
[default  is  ^7  seconds] ; 

5)  POINTFL6  >  initial  status  for  star  field  matching  file 
(FALSE  <■  no  star  field  matching  file); 

6)  PATHFLG  -  initial  status  for  CAS  file  (FALSE  «  no  CAS  file). 

Frooessisg: 

tiNiTt  (Initialization): 

Set  PTSPEC(l)  «  PTSPEC(2)  «  0  [Time  specifications  for  star  field 
pointing  data.] 

Set  ATSPEC(l)  =  ATSPEC(2)  =  0  [Time  specifications  for  CAS  data.] 

Read  input  specifications  [See  "User  Processing  Options"  above.] 

and  convert  start  and  end  specifications  to  VIPER  Julian  day  and 

fraction. 

iVRXNPl: 

Read  VRDATA  items:  VDATE  (date;  YYMMDD) ,  VTIME  (time;  SSSSS.FFF), 
SHPOS  (shutter  position),  FILPOS  (filter  position),  VTENp 
(temperature),  NBUF  (buffer  size),  VRRAD(N) ,  N  =  I,  NBUF 
(radiance) . 

If  end-of-file  then  Go  to  LASTBUF. 

Merge  VDATE,  VTIME  into  day  and  fraction  form:  VTSPEC  = 


46 


Appendix  F  -  SPATIAL  Program's  Processing  Description 


NDAyS({VDATE))  VTIME(in  seconds) /86400  (double  precision). 

If  VTSPEC  <  start  time  then  Go  to  VRINP.  [could  lose  remainder 
of  up  to  1  second  of  data  in  block  after  start  time.] 

If  VTSPEC  >  end  time  then  Go  to  LASTBUF. 

Begin  loop  VRSAMPL,  for  NSAMPL  «  1  to  NBUF: 

Initialize  visual  radiometer  pointing  position,  assuming  no 
pointing  Information  is  available  (can  revise  this  pointing 
later) : 

SPVRRA  «  0; 

SPVRDEC  -  100. 

SRCFLG  *  0  (pointing  source  flag  initialization). 

(STARIMEl: 

If  POINTFLG  -  FALSE  then  Go  to  PATHINP. 

If  PTSPEC(l)  <-  VTSPEC  <  PTSPEC(2)  then  Go  to  POINTING. 

(For  case  of  VTSPEC  outside  PTSPEC  bounds:) 

Read  MATCHING.DAT  items,  searching  for  matching  or  bracketing 
time  for  VTSPEC: 

PDATE(I)  (date:*MMDD) ,  PTIME(I)  (time:  HHMNSS.FF),  VRRA(I) 

(VR  right  ascension,  degrees),  VRDEC(I)  (VR  declination, 
degrees),  TVFLG(I)  (camera),  I  >  1,  2. 

If  end-of-file  for  I  -  1,  then 
Set  POINTFLG  -  FALSE; 

SET  PTSPEC(l)  •  PTSPEC(2)  »  0; 

Go  to  PATHINP; 

Endif 

If  end-of-file  for  I  «  2,  then 

Set  PDATE(2)  -  PDATE(l),  PTIME(2)  -  PTIME(l),  VRRA(2)  » 
VRRA(l),  VRDEC(2)  »  VRDEC(l); 

Set  POINTFLG  »  FALSE; 

Endif 

Merge  PDATE,  PTINE  into  day  and  fraction  form,  using  year 

supplied  in  input  specifications  [This  may  be  performed  while 
searching  for  matching/bracketing  values.]:  PTSPEC(I)  » 
NDAyS({PDATE(I) ))  -t-  PTINE(I)  (in  seconds) /86400  (double 
precision);  should  have  PTSPEC(l)  <-  VTSPEC  <  PTSPEC(2)  [if 
VTSPEC  >  PTSPEC(2)  due  to  EOF,  will  ignore  values  later]. 

(POINTING): 

[Check  that  successive  star  field  references  are  not  so  far 
apart  as  to  be  invalid  references.] 

If  PTSPEC(l)  <  (VTSPEC  -  TIMEINT/86400)  or  PTSPEC(2)  >  (VTSPEC 
+  TIMEINT/86400)  then  Go  to  PATHINP; 

If  PTSPEC(2)  <  VTSPEC  then  Go  to  PATHINP; 

If  VTSPEC  =  PTSPEC(l)  Then 

SPVRRA  =  VRRA(l);  (Visual  Radiometer  right  ascension) 

SPVRDEC  =  VRDEC(l);  (Visual  Radiometer  declination) 

If  TVFLG(l)  =  0  [Pulnix]  Then 
Set  SRCFLG  =  2; 

Else  If  TVFLG(l)  =  1  [Xybion]  Then 
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Set  SRCFL6  -  1; 

Endif 

Else 

Interpolate  between  (VRRA(l) ,VRDEC(1) )  and  (VRRA(2) ,VRDEC(2) ) 
according  to  time  intervals  defined  by  Pl’SPEC*  VTSPEC  to 
obtain  SPVRRA  and  SPVRDEC; 

If  TVFLG(l)  -  TVFLG(2)  «  0  (Pulnix)  Then 
SRCFLG  *  2  * 

Else  If  TVFLG(l)'-  TVFL6(2)  -  1  (Xybion]  Then 
Set  SRCFLG  -  1; 

Else 

Set  SRCFLG  ■  3;  (both  cameras  used) 

Endif 

Endif 

UBAXUim: 

If  PATHFLG  >  FALSE  then  Go  to  NOPATH. 

If  ATSPEC(l)  <-  VTSPEC  <  ATSPBC(2)  then  Go  to  PATH. 

(For  case  of  VTSPEC  outside  ATSPEC  bounds.] 

Read  CAS  it«tms«  searching  for  matching  or  bracketing  time  for 
VTSPEC • 

ADATE(^i)  j  -  i,  3  (date:  year,  month,  day),  ATINE(J,X)  J  » 
1,  3  (vime:  hour,  minute,  second),  AEUL(J,1)  J  •  1,  3 
(Euler  angles  for  shuttle  orientation),  APOS(J,l)  J  •  l,  2 
(right  aseension/declination  angles  for  shuttle  position), 
ARAD(I)  (radial  position  for  shuttle),  ORBXNC(X) 
(inclination  of  shuttle  orbit),  ORBNODE(X)  (right  ascension 
of  shuttle  orbit  ascending  node),  ASUN(J,X)  J  »  1,  2  (right 
ascension/declination  angles  for  sun),  ANOON(J,X)  J  «  1,  2 
(right  ascension/declination  angles  for  moon),  I  »  1,  2. 

Xf  end-of-file  for  X  »  l,  then 
Set  PATHFLG  -  FALSE; 

SET  ATSPEC(l)  »  ATSPEC (2)  »  0; 

Go  to  NOPATH; 

Endif 

Xf  end-of-file  for  X  =  2,  th«m 

Set  A0ATE(J,2)  *  ADATE(J,1],  ATXHE(J,2)  »  ATXME(J,I). 
AEUL(J,2)  »  AEUL(J,1),  A7>OS(J,2)  »  AP0S(J,1),  ARAD(2)  = 
ARAD(l),  ASUN(J,2)  =>  ASry(J,l),  ANOON(J,2)  »  AHOON(J,l); 

Set  PATHFLG  »  FALSE; 

Endif 

Merge  ADATE,  ATIME  into  day  and  fraction  form  (this  may  be 
performed  while  searching  for  matching/bracketing  values): 
ATSPEC(X)  =  NDAYS((ADATE(J,I) ))  +  ATIME( 1, I)/24  + 
ATXME(2,X)/1440  ^  ATXME(3, f )/86400  (double  precision);  should 
have  ATSPEC (1)  <=  VTSPEC  <  ATSPEC(2)  [if  VTSPEC  >  ATSPEC(2) 
due  to  EOF,  will  ignore  values  later]. 

lEMm: 

[Check  that  successive  CAS  references  are  not  so  far  apart  as 
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to  be  invalid  references.] 

If  ATSPEC(l)  <  (VTSPEC  -  TIMEINT/86400)  or  ATSPEC(2)  >  (VTSPEC 
+  TIMEINT/86400)  then  Go  to  NOPATH; 

If  ATSPEC(2)  <  VTSPEC  then  Go  to  NOPATH  (beyond  end  of  CAS 
data ) ; 

If  SRCPLG  .ne.  0  then  Go  to  CELEST2  (can  s)cip  pointing 
determination  from  CAS  if  available  from  star  fields) . 

If  VRZEN  ■  200  then  Go  to  CELEST2  (can  skip  pointing 
determination  from  CAS  if  no  body  coordinate  pointing 
reference  is  available) . 

If  VTSPEC  -  ATSPEC(l)  Then 

Set  SINGLE  «  TRUE  (to  indicate  that  only  first  transformation 
is  needed  to  specify  visual  radiometer  pointing) . 

Use  ORBINC(l),  ORBNODE(l),  AP0S(J,1),  J  -  1,  2,  ARAD(l)  to 
compute  transformation  matrix  TUVWl  for  transformation  from 
shuttle  UVW  coordinates  to  celestial  coordinates  (Procedure 
UVWCEL) . 

Use  Euler  angles  AEUL(J,1),  J  ■  1,  3  to  compute 

transformation  matrix  TXYZl  for  transformation  from  shuttle 
body  coordinates  to  UVW  coordinates  (Procedure  XYZUVW) . 
Compute  pointing  direction  for  visual  radiometer  in  celestial 
coordinates,  for  time  ATSPEC(l)  (Procedure  PATHVR) . 

( CELESTl > ; 

Else 

Set  SINGLE  «  FALSE  (Indicating  that  both  transformations  are 
needed  to  specify  visual  radiometer  pointing.] 

Use  ORBINC(l),  ORBNODE(l),  AP0S(J,1),  J  «  1,  2,  to  compute 
transformation  matrix  TUVWl  for  transformation  from  shuttle 
UVW  coordinates  to  celestial  coordinates  (Procedure 
UVWCEL) . 

(Set  a  flag  to  determine  whether  recalculation  of  the 
transformation  matrix  TXYZl  is  required,  based  on  new  Euler 
angle  inputs.  Also,  test  the  time  to  determine  if 
possibility  exists  of  using  an  earlier  TXYZZ  for  this 
transformation. ] 

Use  Euler  angles  AEUL(J,1),  J  »  1,  3  to  compute 

transformation  matrix  TXYZl  for  transformation  from  shuttle 
body  coordinates  to  UVW  coordinates  (Procedure  XYZUVW). 
Compute  pointing  direction  for  visual  radiometer  in  celestial 
coordinates,  for  time  ATSPEC(l)  (Procedure  PATHVR). 

Use  ORBINC(2),  ORBNODE(2).  AP0S(J,2),  J  »  1,  2,  to  compute 
transformation  matrix  TUVW2  for  transformation  from  shuttle 
UVW  coordinates  to  celestial  coordinates  (Procedure 
UVWCEL) . 

(Set  a  flag  to  determine  whether  recalculation  of  the 
trans:  format  ion  matrix  TXYZ2  is  required,  based  on  new  Euler 
angle  inputs.] 

Use  Euler  angles  AEUL(J,2),  J  »  i,  3  to  compute 

transformation  matrix  TXYZ2  for  transformation  from  shuttle 
body  coordinates  to  UVW  coordinates  (Procedure  XYZUVW). 
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Compute  pointing  direction  for  visual  radiometer  in  celestial 
coordinates,  for  time  ATSPEC(2)  (Procedure  PATHVR) . 
Interpolate  pointing  directions  for  visual  radiometer  at 
times  ATSPEC(l)  and  ATSPEC(2)  to  obtain  visual  radiometer 
pointing  at  time  VTSPEC. 

Endif 

ismsnii 

[Store  celestial  object  and  shuttle  reference  positions,  after 
interpolating  between  reference  times  ATSPEC(l)  and  ATSPEC(2), 
if  necessary.] 

(SUNRA,  SUNDEC)  -  Interp(ASUN,  ATSPEC,  VTSPEC,  SINGLE) 

(LUNRA,  LUNDEC)  >  lnterp(AMOON,  ATSPEC,  VTSPEC,  SINGLE) 

SRCFL6  *  8RCFL6  4  (retain  camera  attribute;  indicate  CAS  data 
is  present) . 

Go  to  STOREREC. 

( NOPATH >t  (set  fill-in  values) 

ORBALT  •  -1; 

ORBRA  -  SUNRA  -  LUNRA  »  0; 

ORBDEC  -  SUNDEC  «  LUNDEC  »  100. 

ISIQREREC) : 

Store  record  in  output  buffer. 

Increment  NRECS. 

If  NRECS  -  LIMBUF  then 
Write  output  data  block. 

Set  NRECS  «  0 
Endif 

Increment  time  for  visual  radiometer  sample:  VTSPEC  =  VTSPEC  -f 
0.01/86400. 

Next  NSANPL. 

End  VRSANPL  loop. 

Stop. 


( LASTBUF ) ! 

If  NRECS  .ne.  0  then  Write  output  data  block  (NRECS  <=  LIMBUF). 
Stop. 


(UVWCELK 

A  s  APOS(l,I)  (shuttle  position  right  ascension); 

0  »  APOS(2,I)  (shuttle  position  declination); 

K  B  ORBINC(I)  (shuttle  orbit  inclination)  (REAL); 

N  »  ORBNODE(I)  (shuttle  orbit  ascending  node)  (REAL) ; 


0  =  (Ul,  U2,  U3)  =  (cos  A  cos  D,  sin  A  cos  D,  sin  D) ; 

W  =  (Wl,  W2,  W3)  =  (sin  K  sin  N,  -sin  K  cos  N,  cos  K) ; 

V  =  (VI,  V2,  V3)  =  cross-prod(W,  U)  (vector  product); 
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ui  VI  wi 

TUVWI^m  V2  W2 
'■  [U3  V3  W3_ 

f.  Return  TUVWI . 

immit 

A  >  AEUL(1|I)  (shuttle  yaw  Euler  angle); 

B  »  AEUL(2,I)  (shuttle  pitch  Euler  angle); 

P  »  AEUL(3,I)  (shuttle  roll  Euler  angle); 

10  0 
TX-  0  COS  P  SIN  P 
,0  -SIN  P  COS  P 

COS  B  0  -SIN  B 
7Y-  0  1  0 

SIN  BO  COS  B 

COS  A  SIN  A  0 
rZ-  -SIN  A  COS  A  0 
0  0  1 

TXYZI  -  TX  *  TY  *  TZ 
Return  TXYZI. 

(PATHVR>; 

A  »  VRAZ  (VR  azimuth,  in  shuttle  coordinates) ; 

Z  »  VRZEN  (VR  zenith,  in  shuttle  coordinates) ;  [nominally  180 
degrees] 

L  =  (cos  A  cos  Z,  sin  A  cos  Z,  sin  Z) ; 

[The  above  three  steps  are  only  required  for  the  initial  call  for 
a  given  data  processing  run,  because  VRZEN  and  VRAZ  remain 
f ixed . ] 

1  SPVRI  =  TUVWI  *  TXYZI  *  L; 

[Decompose  SPVRI  to  right  ascension  and  declination  as  follows: 
SPVRRAI  =  atan2(  SPVRI2,  SPVRIl)  (degrees); 

,  SPVRDECI  =  acos  (  SPVRI3)  (degrees).] 

Return. 
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